Power Spectra in Global Defect theories of Cosmic Structure Formation 
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An efficient technique for computing perturbation power spectra in field ordering theories of cosmic 
structure formation is introduced, enabling computations to be carried out with unprecedented 
precision. Large scale simulations are used to measure unequal time correlators of the source stress 
energy, taking advantage of scaling during matter and radiation domination, and causality, to make 
optimal use of the available dynamic range. The correlators are then re-expressed in terms of a 
sum of eigenvector products, a representation which we argue is optimal, enabling the computation 
of the final power spectra to be performed at high accuracy. Microwave anisotropy and matter 
perturbation power spectra for global strings, monopoles, textures and non-topological textures are 
presented and compared with recent observations. 



I. INTRODUCTION 

Over the next few years high resolution maps of the 
cosmic microwave sky will become available. These maps 
will allow competing theories of cosmic structure forma- 
tion to be tested with exquisite precision. The primary 
quantities of interest for comparing theories to observa- 
tions are the power spectra of fluctuations, both for the 
microwave sky temperature and the mass density. The 
predictions of simple inflationary theories have largely 
been worked out, typically to one per cent accuracy, and 
detailed comparisons with data are taking place. The 
state of the main rival set of theories, based on sym- 
metry breaking and phase ordering, has been less rosy. 
These theories involve a causal source comprising the or- 
dering fields and/or defects, which continually perturbs 
the universe on ever larger scales. In the inflationary the- 
ories, linear perturbation evolution is all that is needed: 
but for the defect theories full linear response theory is 
required. The defect sources are in general 'decoherent' 
0, providing additional computational difficulty. This 
Letter presents a general solution method for solving the 
linear response problem in such models. 



II. STIFF SOURCES: MEASURING UNEQUAL 
TIME CORRELATORS 

Accurate codes have been developed for field evolu- 
tion in different symmetry breaking theories. But a sin- 
gle simulation cannot simultaneously resolve all scales 
of observational relevance. Thus information gathered 
from simulations must be combined in some way. It has 
been clear for some time that the ideal quantity which 
a) uses all the information present in a simulation, b) 
incorporates the powerful properties of scaling evolution 
and causality, and c) preserves all the information needed 
to compute power spectra is the unequal time correlator 
(UETC) of the defect source stress energy tensor Q^v'- 



(e^,(k, r)epA(-k, r')) = Cp.,pA(A:, r, r') 



(1) 



where t, r' denote conformal time, and k comoving 
wavenumber. Because the perturbed Einstein-matter 
equations are linear, all perturbations are determined 
in terms of the source via appropriate Greens functions. 
Thus in principle all quadratic estimators of the density 
perturbations are determined by (Q). 

The unequal time correlators are highly constrained by 
causality, scaling and stress energy conservation. Causal- 
ity means that the real space correlators of the fluctuat- 
ing part of must be zero for r > r + r' Scal- 
ing ||] dictates that in the pure matter or radiation eras 
C^l,y,p\ oc 4>q/ {tt')2 Cf,^^p\{kT, kr'), where i^o is the sym- 
metry breaking scale and c is a scaling function. Finally, 
energy and momentum conservation for the stiff source 
(see e.g. [^) provide two linear constraints on the four 
scalar components of the source. Any pair determines 
the other two up to possible integration constants. In 
our work we have found the best pair to be the energy 
density Oqo and the anisotropic stress Q^; the energy 
and momentum conservation equations give good behav- 
ior for all components on both superhorizon and subhori- 
zon scales. This choice is also favoured by the fact that 
000 and , along with the vector and tensor compo- 
nents, and O"'" fix all superhorizon perturbations in 
the most direct manner. However, we have also checked 
that other choices give consistent results. 

As mentioned, our method uses scaling and causality to 
extend the dynamic range of the numerical simulations. 
In the simulations, the fields start from uncorrelated, but 
non-scaling initial conditions, and evolve towards scaling 
behaviour. We evolve the fields to some final time r 
when the system is well into the scaling regime. At this 
time we compute the stress energy tensor Q^i,, which 
we Fast Fourier Transform and decompose into the vari- 
ables 9oo, 0'^, 0^, and 0"^. We then repeat the sim- 
ulation with identical initial conditions, computing the 
same quantities at several times t' < t. This procedure 
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ensures that the equal time correlations, which dominate 
the perturbation production, are measured when the sys- 
tem is closest to scaling. We store the isotropic averages 
of the correlators as functions of the wave vector magni- 
tude k, e.g. Coo,s{kT,T'/T) =< eoo(k,T)e'5(-k,r') >. 
Statistical and sampling errors are small at large k be- 
cause many modes contribute, but at small k the sam- 
pling is sparse and the noise larger. In particular the 
value at fc = represents the mean defect density and 
should be discarded. In order to increase the resolution 
for small fcr, we compute the correlators in real space 
and project them onto the radial separation r. The fact 
that the correlators must vanish for r > t + t' allows 
us to solve for an additive constant, in effect determin- 
ing the k — Q value. The projected real space correla- 
tor is then transformed back to Fourier space using the 
continuum Fourier transform over the compact support 

C(k) = 47r /o^^^ r^dr C{r)sm{kr)/kr. For small values of 
k we use the latter equation, while for large values we use 
the direct Fourier space computation. For intermediate 
kr ^ b the two match well, but this procedure allows us 
to extend our dynamic range in small fcr by an order of 
magnitude. Note that using causality in this way limits 
the final time in the simulation to r < j where L is the 
box size; typically we evolve in a 400^ box to r = ^ = 50 
so that there are still many independent horizon volumes 
in the box. 

III. EIGENVECTOR PRODUCT 
REPRESENTATION 

We have explained how one can devote the full numeri- 
cal power at hand to compute the UETCs alone, without 
wasting computational effort or storage on linear pertur- 
bation theory. Next we shall show how very fast Einstein- 
Boltzmann solvers recently developed ||^ can be used to 
convert the UETCs into cosmological power spectra with 
relatively small numerical effort and very high precision. 

This is done by representing the UETCs as a sum of 
eigenvector products [|j . The idea is to regard the stress 
energy correlators (|l|) as 'symmetric matrices' with in- 
dices iJ,iy,T and pXr'. In practise, to compute the scalar 
perturbations we need the auto- and cross-correlators of 
two components (for example Qqo and 6'^), and for the 
vector and tensor perturbations we need the two auto- 
correlators of 9^ and 6"^. Regarded as matrices, the 
correlators involved are symmetric and positive definite 
(expectation values of squares) and so the eigenvalues 
are all real and positive. Matrix index summation is re- 
placed with an integral / dTw{T) with w{t) some chosen 
weighting function. Often the choice of weighting is nat- 
urally dictated by scaling and dimensional analysis, but 
in any case the results were checked to be independent 
of it. For sensible choices of w{t) the trace of the corre- 



lators is finite: it follows that there is a convergent series 
of positive eigenvalues which may be used to index the 
eigenvector sum. The correlators can then be expressed 
as an infinite sum 

C(fc,T,T')=^AV(fc,rK(fc,r'), (2) 

i 

where 

J dT'C{k,T,T'yik,T')w{T') = \'v\k,T). (3) 

The indices labelling the components of the stress tensor 
are implicit. We have found that including the largest 15 
eigenvectors typically reproduces the unequal time corre- 
lators to better than ten per cent precision: the effect of 
including more than 15 on the final perturbation power 
spectra is negligible at the few per cent level. 

We believe the properties of this representation make 
it ideal for the purpose of computing cosmological per- 
turbations. Namely a) the representation automatically 
minimises the 'least squares' fit, for a given number 
of eigenvectors b) the eigenvectors individually conserve 
stress energy c) the eigenvectors individually vanish as 
T goes to zero since the correlators C{k,T,T') vanish for 
T -^t' , making specification of the initial conditions triv- 
ial, and d) the contribution of successive eigenvectors to 
the perturbations converges quickly, as a result both of 
the decrease in eigenvalue and the increasingly oscilla- 
tory nature of the eigenvectors. In particular the inco- 
herent superhorizon (fcr < 1) portion of the correlators 
is represented by an infinite sum of ever more oscillatory 
eigenfunctions, which have increasingly little effect on the 
perturbations. For the computation we discretise the un- 
equal time correlators sampled in 256 equal logarithmic 
spacings in the interval e~^ "' < fcr < e^'^. The individual 
eigenvectors are each fed into a full Boltzmann code [||, 
and the total perturbation power spectra are then given 
by the sum of those for individual eigenvectors, weighted 
by their eigenvalue. 

During the pure matter and radiation epochs, the pro- 
cedure is simplified because the correlators scale, and 
so are represented for all /c by a single set of eigen- 
vectors, functions of kr. We incorporate the matter- 
radiation transition by repeating the computation of un- 
equal time correlators for several different values (typ- 
ically 20) of t/t*, where r* is the conformal time at 
matter-radiation equality ^^.(t*) = U,m{T*). After the 
simulations are completed, we collect the results in a 
matrix of correlators, H(fcT, fcr', fcr*), which are then di- 
agonalised to produce a set of eigenvectors for each fcr* 
considered. These eigenvectors smoothly interpolate be- 
tween those for the radiation era and those for the mat- 
ter era, so the Boltzmann code can use a simple spline 
interpolation between them. The integration solves the 
full linearised relativistic Einstein equations tracking the 
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FIG. 1. Angular power spectrum of anisotropies generated 
by the scalar component of the source stress energy for global 
strings. The upper curve shows the total spectrum, the lower 
ones contributions from individual eigenvectors. This Fig- 
ure illustrates decoherence: each eigenvector individually pro- 
duces an oscillatory Ci spectrum, but these oscillations all 
cancel in the sum. 

photons 5r, baryons 5b, cold dark matter 5c, neutrinos 
5u, and their relative couplings. The evolution includes 
the full matter-radiation transition, finite recombination 
rate, and other effects. The Boltzmann calculations are 
accurate to about 1%. 

IV. CHECKS 

Many checks have been performed on this procedure, 
which we briefly summarise here Q| . When the box size of 
texture simulations is reduced from 400'^ to 256"^, the C; 
spectrum changes by less than 3%. To check consistency 
with energy conservation we used a different pair of scalar 
variables 8oo + 0, 8 + 28'^, which changed the results by 
5%, in 400'^ boxes, but the latter pair gave much slower 
convergence with box size. Weighting the diagonalisa- 
tion by an additional factor of t^^^ adds more weight to 
subhorizon scales, and affects the result by up to 10% at 
I = 2, but less than 2% at larger We have compared 
the total, as well as scalar, vector and tensor anisotropies 
to those produced by the direct line-of-sight integration 
code and they individually agree to 10-20%. Since the 
current method includes additional contributions at the 
last scattering surface, explicitly uses scaling, employs 
far larger box sizes, and more accurate integrations, we 
conclude the results are consistent with the new results 
being much more accurate. All these checks indicate that 
these new results are reliable to better than 10% [Q. 



10 100 1000 10 100 1000 



' total ' ' " ' . 

scalar / \ strings , 

vector / \ 

tensor/ \ 

^\ \ 

-^^-.y - 


1 1 1 

^ — V monopoles . 

1 1 Mi^H 


1 1 1 — 1 1 1 — 1 iTrwii| • 

textures 

^ ^ \ 
^ \ 

.V...... 


nt. textures 

1 1 ■ 1 



10 100 1000 10 100 1000 

I I 

FIG. 2. The contributions to the total power from scalar, 
vector and tensor components. 

V. RESULTS 

The resulting angular power spectrum of the temper- 
ature anisotropies for the scalar perturbations sourced 
by global cosmic strings are shown in Figure |. Each 
of the eigenmodes is coherent and shows the expected 
acoustic peaks Decoherence is manifested when we 
sum over eigenvectors. The peaks add out of phase, 
resulting in a smooth angular power spectrum. One 
expects cosmic strings to be the least coherent of the 
defects under consideration. Indeed, the = 6 non- 
topological texture scalar modes do exhibit residual oscil- 
lations even when summed over all eigenmodes as shown 
in Figure |^. The figures show C; spectra computed for 
h = 0.5, rib = 0.05, f2 = 1. The dependence on the 
Hubble constant h and baryon content fib is weak 
The most striking feature of all the models under con- 
sideration is the predominance of vector modes. They 
dominate up to Z 100, at which point they are sup- 
pressed by the horizon size on the surface of last scat- 
ter. It is not hard to see that vector and tensor con- 
tributions should be at least comparable to the scalar 
contribution. By causality the k space correlator is an 
integral over a real space function of compact support, 
and should be analytic in k at small k. By isotropy, we 
can expand any correlator at small kT as Cij ^^i (kr, kr') = 
A5ij5ki -t- B{5if^5ji + 5ii5jk) + 0{k^). The constant A con- 
tributes only to trace correlators: the single constant B 
then determines the anisotropic scalar, vector and tensor 
contributions, in the ratios cs.s ■ cv,v : ct.t = 3:2:4 
(these are accurately verified in our simulations). Thus 
vectors and tensors can be expected to contribute a sig- 
nificant fraction of the temperature anisotropies in field 
ordering theories. (A general analysis of the relative im- 
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FIG. 3. Comparison of defect model predictions to current 
experimental data. All models were COBE normalised at 
I = 10. 



portance of vector and tensor modes will be described 
elsewhere Q.) The large amplitude of vector modes and 
the decoherent sum of eigenmodes leads to a suppres- 
sion of power at / ^ 100 a very different spectrum to 
that expected from adiabatic fluctuations in inflationary 
models. We show a comparison between the predictions 
of the global field defect theories and the current gener- 
ation of CMB experiments in Figure |^. All models are 
normalised to COBE at I — 10. They are all systemat- 
ically lower than the current degree-scale experimental 
points. 

The same calculations directly yield the matter power 
spectrum shown in Figure ^ Normalised to COBE, our 
tests indicate that the results should be reliable to a few 
percent. From the power spectra we derive the nor- 
malization as of the matter fluctuation in 8h~^ Mpc 
spheres. Global strings, monopoles, texture and iV = 6 
non-topological texture give erg = 0.26, 0.25, 0.23, 0.21, 
respectively, for h = 0.5, and scaling approximately as 
h. The field normalization for textures is e = 8tt'^G(I)q ~ 
1.0 X 10"'*, consistent with our previous calculation ||] 
of e = 1.1 X lO"**. These normalizations are a factor of 
5 lower than the generic prediction of n = 1 inflationary 
models where as = 1.2 for h = 0.5. Cluster abundances 
suggest values of cts ~ 0.5 for a flat universe. 

To summarise, the techniques used here enable us 
to convert unequal time correlators into temperature 
anisotropy and matter fluctuation power spectra within 
a few hours on a workstation. For all the defect theo- 
ries, vectors contribute approximately half of the total 
CMB anisotropy on large scales, leading to a suppression 
of acoustic peaks and a low normalization of the matter 
power spectrum cg ~ O.25/150. Current observations of 
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FIG. 4. Matter power spectra computed from the Boltz- 
mann code summed over the eigenmodes. The upper curve 
shows the standard cold dark matter (sCDM) power spec- 
trum. The defects generally have more power on small scales 
than large scales relative to the adiabatic sCDM model. The 
data points show the mass power spectrum as inferred from 
the galaxy distribution fetj. 



CMB anisotropics and galaxy clustering do not favor the 
models under consideration. 
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